Full counting statistics for SU(A^) impurity Anderson model 
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We analyze the full counting statistics of a multiorbital Kondo effect in a quantum dot with the 
SU(A'^) symmetry in the framework of the renormalized perturbation theory. The current prob- 
ability distribution function is calculated for an arbitrary dot-site Coulomb repulsion U in the 
particle-hole symmetric case. The resulting cumulant up to the leading nonlinear term of applied 
bias voltages indicates two types of electron transfer, respectively carrying charge e and 2e, with 
different A''-dependences. The cross correlation between different orbital currents shows exponential 
enhancement with respect to U, which directly addresses formation of the orbital-singlet state. 
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Experimental realization of the Kondo state in meso- 
scopic devices has promoted further research of the 
Kondo effect since it enables one to access many-body 
effects in out-of-equilibrium under finite bias voltagesJ^ 
Recent theoretical studies of nonequilibrium noise for 
Kondo dota^"— have predicted enhancement of shot noise 
with a fractional Fano factor, and have stimulated sub- 
sequent experimental studieSfiSrJ^ 

Deeper understanding of the nonequilibrium Kondo ef- 
fect is brought in by the current distribution function, 
which includes higher-order cumulants beyond the first 
and second cumulant (average current and noise power). 
It, however, is still challenging to calculate the cumulant 
generating function (CGF) in the nonequilibrium Kondo 
states. Recently, Komnik, Gogolin and Schmidt have de- 
rived current probability distribution for the SU(2) An- 
derson impurity^"— from the general formulation of the 
full counting statistics (FCS)] ^^^^ They have clarified 
that the nonequilibrium backscattering current is com- 
posed of two types of electron transfer due to a single 
quasiparticle and a pair of quasiparticles, carrying charge 
e and 2e respectively. 

In the present work, we extend the FCS approach to 
multiorbital Kondo dots, which have been experimentally 
investigated in vertical dots,^^ carbon nanotubes,^^ and 
double dotsj^ We consider an SU(7V) impurity Anderson 
model as a prototype model for examining multiorbital 
effects. By employing the renormalized perturbation the- 
ory (RPT)^^'^^, we calculate the zero-temperature CGF 
for the entire strength of the dot-site Coulomb repulsion. 
The RPT is based on a general idea of renormalization 
in quantum field theory, and is consistent with several 
known results, indicating its validity up to terms with 
third order of applied bias voltage. ^-i^ Our calculation 
provides direct information on an orbital singlet, i.e., a 
correlated electronic state involving different orbitals. 

Model — Let us consider a single quantum dot system 
described by the SU(iV) impurity Anderson model = 
Ho + "^T + "^c/ with 



(2) 
(3) 



where d^^ annihilates an electron in the dot level with 
orbital m — 1,2, ■ ■ ■ , N , Cf.^^ annihilates a conduction 
electron with moment k and orbital m in lead a = L, R, 
and U is the dot-site Coulomb repulsion. Here, orbital 
includes spin, and thus N is even. The intrinsic level 
width of the dot levels owing to tunnel coupling Va, is 
given by r = J^a ''^Pc\va\'^ and the density of state of the 
conduction electrons pc- For simplicity, the symmetric 
lead-dot coupling Vj^ — and the particle-hole sym- 
metry — ~[N — l)C^/2 are assumed. The chemical 
potentials Ml/_r — ±^2, satisfying — /i^ = V{> 0), 
are measured relative to the Fermi level which is defined 
at zero voltage V = 0. We take a unit, Ti = ks = e ~ 1, 
throughout this paper. 

Full counting statistics — The probability distribution 
P{q) of the transferred charge q = [qi,q2, - ■ ■ ,qj^) with 
orbital subscript across the dot during a time interval T 
provides current correlation function of all orders. In this 
paper, we define the transferred charge operator as q^^ = 
"Lm(-7'/2) - n^„(r/2), where n^^^it) is the electron 
number in the left lead. In order to discuss the correlation 
functions systematically, we calculate the CGF Inx (A) = 
In^g e'^^''^P{q) in the Keldysh formalism,^- 



lnx(A) = ln(rc^^), 



(4) 



where = Tcexp{-i [H^ (i) (i)] } is the 
time evolution operator for an extended Hamiltonian 
Ti^ ~ Hq + + Hjj, C is the Keldysh contour along 
[t : -r/2 ^ +T /2 -T /2], Tc is the contour ordering 



operator, and A = {X^, A2, ■ 
Here, is given by 



A^) is the counting field. 
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with the contour dependent counting-field defined by 



^0 = E ^k^^ncLAr.. + E ^d^rfld^ , (1) Kit) = A ^ ±\„ for the forward and backward 



paths labeled, respectively, by 



and 



2 



In order to calculate the CGF in Eq. (|4]), we make 
use of a procedure suggested by Komnik and Gogoliuj^ 
which is outlined below. First, a more general function 
x(A_, A_,_) with A-|- = (A]^^, Aa^, • • • , Ajyzp) is introduced. 
It is basically given by Eq. (|4]) but A^^ is treated formally 
as an independent variable assigned for each contour. For 
the long time limit T —> oo where the switching effects 
are negligible, the general CGF is proportional to T, 



Inx 



(6) 



with the adiabatic potential U{\_,Xj^). Once the adi- 
abatic potential is computed, the statistics is recovered 
from Inx(A) = —iT'U^X, — A). Performing the derivative 
of Eqs. (ID) and © with respect to A„j_ , we obtain 



-iY(A_,A+) ^ lim 



d 



dX„i~ T-i-oo \ dXr. 

where we use a notation 



-H^(O-)) ,(7) 



A .A_, 



(r^5^A(i))/x(A_,A+). 



(8) 



Equation ^ presents an expectation for Hamiltonian H\ 
and the Wick's theorem is applicable. The right hand 
side of Eq. ([7]) can be expressed in terms of the Green's 
function. 



dXr 



-W(A_,A+) 



duj 
2^ 



iA^/2 — h 



5Lfe™MGiirM ,(9) 



with A„ = A,„_ - A„^+. .gfc„ + (w) = i2TrS{uj - 

^kam)fJ^) and = -i27td{uj - £fe„„)[l - 

are the lesser and greater parts of the Green's func- 
tion for electrons in lead a with the Fermi distri- 
bution function fa{uj) = [e^""^")/^ -|- respec- 
tively. For the long time limit T oo, the dot 
Green's function is defined as {G'^„j(a;)} , = —i J d{t — 

f)eMt-t')(^TcdMdlK,))^_^^^. Here% and i^' are 
the labels for the two Keldysh contours. 

Renormalized perturbation theory — The three basic 
parameters that specify the impurity Anderson model 
are edm, F and U . The low-energy properties can 
be characterized by the quasiparticles with the renor- 
malized dot-level 'cdm = z [edm + '^^dmi^)]^ renormalized 
level width T = zT, and the renormalized interaction 



U 



2r(4) 



(0,0; 0,0) (to ^ m'), where T^dmi^) is 
the self-energy of the retarded Green's function for the 



dot state: C^dmi^) 



dm 



[1 — 9EJ^(w)/9w|;j=o] ^ is the wave function renormal- 

ization factor, and T\J^,{u}i,uj2',^z,'-^i) is the local full 
four-point vertex function for the scattering of the elec- 
trons with the orbital to and to'i^i^ Note that these 
parameters are defined at equilibrium V = Q and A = 0. 



The replacement of the bare parameters with the renor- 
malized ones gives leading terms of the Hamiltonian cor- 
responding to the low-energy fixed point of the Anderson 
model in Wilson's theory,^'' The perturbation theory in 
powers of U can be reorganized as an expansion with re- 
spect to the renormalized interaction J7, taking the free 
quasiparticle Green's function ^^^^^(a;) = (w — edm + ^F)^^ 
as the zero-order propagator. In addition, three counter 
terms are introduced in order to prevent overcounting. 
This procedure has enabled one to calculate the exact 
form of the Green's function in the absence of the count- 
ing fields Gdmi'^)\\=o at low energies up to terms of 
order w^, V^, and ^^24 'pj-^gj-g ^re explicit relations 
between the renormalized parameters and the enhance- 
ment factor of susceptibilities zXd = l + f^P^m(0)j ^Xc d ~ 
\ — {N — l)C^Prf,„(0), and the Friedel's sum rule 7rn^^ — 
cot~^ {^dm/^^ with the electron occupation in orbital to 
of the dot = (dl^d^^), and the renormalized density 
of state Pdm{oj) = (f /7r)/[(w - e^^)^ + f\ In partic- 
ular, for the particle-hole symmetric case (f^^„ = 1/2), 
the renormalized parameters can be expressed in a sim- 



ple form as, [//(ttF) 
Wilson ratio R 



R 



N 



with the 
r can be 



1 and edm = 
^ - N/[{N - 1) + x. Jxd 
considered as the Kondo temperature as, Tk = ttT/4. 
In this paper, we evaluate the renormalized parameters 
by the Bethe ansatz exact solution (BAE) ^^'^^ and the 
numerical renormalization group (NRG) calculation2Zi2£. 

Let us now apply the RPT to calculation of the dot 
Green's function G^,^(a;) for the extended Hamiltonian 
The dot Green's function is given by 

GL(^) = zG'dM = z IgLi^f' - SL(^)1 ~\m 



where the zero-order part is given byi* 

9^--i^) =L + tf ih - 1/2) + if ifR - 1/2) 



'V 



9dnM^) 

~9di-{-) 
9diH^) 

with 



'V 



-uj + lY ifL - 1/2) + zF ifR - 1/2) 



"Dm, 
(11) 



p„(c.) = c.2 + r2 + r2 



-1 (1-/l)/^ 



+ (e'^"/2-l) (l-/^)/4l2) 



The remainder part of the renormalized self-energy can 
be readily calculated in the second order perturbation in 
[/, at T = up to uj'^^ujV and V'^, as 



^dmi^) 




-B*^{-oj,-V) An(a;,F) 



(1?) 



3 



with 



(N - 1) [a (w, 3V/2) + 3 a (w, V/2) 

+ 3a{oj, -V/2) + a {u, -3V/2)] , (14) 

\^e-''^n.'/^(uj,3V/2) 



E 



0/2 



b{io,V/2) 
b{uj,~V/2) 



[2 + e<^" 

+ e'(^-+^"')/26 (w, -3V^/2)} (15) 

Here, a{uj,x) — — i(a; — a;)^sgii(— cj + a;), and, b{uj,x) — 
{uj — x)^9{—uj + x). For A = 0, the exact ex- 
pression of usual self-energy up to uj^,ujV and V"^ is 
reproduced,^i2i Substituting Eqs. ^ and ^ into 
Eq. (jlOp . we readily obtain an expression for the Green's 
function up to uj'^,ijjV and V^. 

Results and discussion — Substituting the obtained 
Green's function into Eq. ([9]) and integrating them re- 
spect with A^_, the CGF is derived at T = up to y'^, 
as Inx(A) = Jq + -^1 



J-2 with 



7- fV/^ 

247r r2 ^ 



-i\r, 



1 



(17) 



J"2 = 



- 1 



(18) 



takes sum of all combination of m and m' 
is the GGF of the zero-order part 



Here, '}2[m^m') 
without m — m 
of the RPT, is the sm^fe-quasiparticle backscattering 
process carrying charge e, and is the <wo-quasiparticle 
backscattering process where the two quasiparticles with 
different orbitals make a singlet state, carrying charge 
2e. In particular, and J^2 represent the reflection 
of quasiparticles by the residual interaction U. This 
CGF for arbitrary U corresponds to the SU(A^) exten- 
sion of the hypothesis presented by Gogolin, Komnik and 
Schmidt. i^J^ii^ 

We now consider the cumulant for full current C„ — 
(— i)"^^ Inx(A), which is derived from the CGF with 
A„j — A for all m, as 

C„ = T [luhn + (-l)"(no + Pbl + 2"Pfc2)] , (19) 

where = NV/{2t:) is the linear-response current and 
5nn' is the Kronecker's delta. We can obtain Pto = 
and ni = m^}MM^Y^ from and F,, 

247r p2 "-^ 247r p2 ^> -1 ' 

respectively. These express the probability of the single- 
quasiparticle backscattering processes carrying charge 
e, per time. Similarly, that of the two-quasiparticle 

backscattering process Pb2 — ^^^~Vi-"~"^^ ^ is ob- 
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FIG. 1. (Color online) (a) The FFIR C^/C^ for n = 2, 
(b) n = 3, and (c) n = 4, as a function of the Coulomb 
repulsion U . (d) The cross cumulant between fluctuation of 
transmitted charge with orbital m and m! m) as a function 
of the Coulomb repulsion U . The blue solid line and the red 
broken line present the SU(2) (A^ = 2) case and the SU(4) 
(A'' = 4) case, respectively. 



with the different orbitals form a singlet pair and carry 
charge 2e, which causes the factor 2" in Eq. 

In order to extract universal properties from the cu- 
mulant (I19p . we consider the Fano-factor inspired ratio 
(FFIR)^'^ for n > 2 and an arbitrary U, normalized in 
the form 



C„ _ l + (l + 2"+i)(iV-l)(fl^-l)^ 
CP l + 5(7V-l)(i?^-l)2 



(20) 



r2 



tained from J-2. In this process, the two quasiparticles 



Here, = (— 1)"/^?" (n > 2) is the Poisson value of 
the CGF and h = In - C^/T = Pto + Ai + 2Pf,2 is 
the backscattering current. Thus, the FFIR can be in- 
terpreted as an average of the cumulants for the three 
backscattering current in this case. Remarkably, the 
FFIR is determined only by the two parameters, the 
Wilson ratio i?^ and degeneracy N. The Fano fac- 
tor which corresponds to the FFIR for n = 2, agrees 
with the previous worksi^^— Evaluating i?^ with the 
BAE for SU(2)(iV = 2)^*29 and with the NRG for 
SU(4)(Ar = i)^L30 we plot the FFIR for noise (n 2), 
skewness {n = 3) and sharpness (n = 4) as a function 
of the Coulomb repulsion U in Figs. [IIa)-[lJc), respec- 
tively. In these figures, with increase of U, the FFIR 
crossovers from the Poisson value to the universal value 
of the SU(2) and SU(4) strong-coupling limit. Note that 
the n dependence, which enters through the factor 2" 
in Eq. (|19p . is caused by the two-quasiparticle process. 
It makes the FCS of the Kondo systems quite different 
from that of the noninteracting system U = 0. In the 
weak coupling limit U ^ 0, the Wilson ratio takes the 
value of Rj^ — 1, and thus the FFIR goes to C^/C^ — s- 1. 
In the strong-coupling limit C/ — > oo, the Wilson ratio 
approaches to the value i?^ N/{N — 1), and the FFIR 



4 



TABLE I. The FFIR for n = 2, 3 and 4 in the strong-couphng 
Umit [/ — !> oo, given in Eq. (|21|) . for several choices of degen- 
eracy A'^. 
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takes a universal form 



C„/C,r^(^ + 2"+i)/(7V + 4). 



(21) 



The explicit values of the FFIR for several N are given in 
the TABLE [J For iV = 2, Eq. (gT]) agrees with the result 
given by Gogolin and Komnik.^"^ In the limit of large 
degeneracy N — > 00, however, the FFIR approaches to 
the Poisson value C„/C^ — >■ 1, even though the Coulomb 
repulsion has been taken first to be [/ — >■ 00. This is 
because the renormalization becomes weaker for larger N 
and the two-quasiparticle process is suppressed Pb2 — ^ 0. 

We next consider cross cumulant which is observed as 
cross correlation between different orbital currents, and it 
may enable one to directly observe the two-quasiparticle 
scattering in experiments. The generic form of the cross 
cumulant between orbital currents is derived from the 
CGF up to terms of as 



{{5q^Jql,))^{-^) 



k+l 



Qk+l 



■Inx(A) 



(22) 



for k,l > I and m ^ m' . The cross cumulant for 
{k,l) = (1,1) is plotted as a function of the Coulomb 
repulsion U in Fig. [T] (d). In the strong-coupling re- 
gion, the cross cumulant is inversely proportional to the 
square of the renormalized level width as {{Sq^Sql^,)) oc 

1/r^, and, thus, increases exponentially with increase 
of the Coulomb repulsion as shown in Fig. [T] (d). The 
positive cross cumulant is a signature of orbital-singlet 



states traveling through the dot. Larger orbital degen- 
eracy makes renormalization weaker and suppresses the 
two-quasiparticle scattering. Therefore, the cross cumu- 
lant becomes the largest in the SU(2) case and the or- 
bital degeneracy suppresses the correlation. We note 
that there is no higher order cross cumulant such as 
{{Sq-in^q'^,5ql^^„)) for m ^ m' ^ m" at low bias voltages 
determined by the terms up to order in the particle- 
hole symmetric case even in the presence of the orbital 
degeneracy N > 2. Naively, it seems that the cross cor- 
relation can be observed in double dots with interdot 
Coulomb repulsion, or dots connected to ferromagnetic 
leads with opposite polarizations. 

Finally, we comment on cumulants of the current for 
orbital m, C™ = (-i)"^lnx(A) = (-1)"/^ with the 
backscattering current /™ = V/h — C™. The cumulant 
C™ always takes the Poisson value jC^^ = 1 in the 
particle-hole symmetric case for the contribution up to 
because there is no scattering of two quasiparticles 
with the same orbital in this case. 

Summary — We have investigated the FCS of a multi- 
orbital Kondo dot described by the particle-hole symmet- 
ric SU(Ar) impurity Anderson model. Using the RPT, we 
derived the CGF for arbitrary Coulomb repulsion up to 
terms of order . The dot-site Coulomb repulsion in- 
duces quasiparticle's orbital-singlet pairs carrying charge 
2e in the backscattering current. This process character- 
izes quantum fluctuations of the current in the correlated 
dot and is particularly manifest in the cross correlation 
between orbital current. It is also found that there is no 
electron entangled state carrying more than three quasi- 
particles in the particle-hole symmetric case in current 
up to order , even in the presence of large orbital de- 
generacy N > 2. 
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